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ABSTRACT 



Context. Recent observational results found a bend in the Tully-Fisher Relation in such a way that low mass systems 

lay below the linear relation described by more massive galaxies. 

Aims. We intend to investigate the origin of the observed features in the stellar and baryonic Tully-Fisher relations and 

analyse the role played by galactic outflows on their determination. 

Methods. Cosmological hydrodynamical simulations which include Supernova feedback were performed in order to follow 

the dynamical evolution of galaxies. 

Results. We found that Supernova feedback is a fundamental process in order to reproduce the observed trends in the 

stellar Tully-Fisher relation. Simulated slow rotating systems tend to have lower stellar masses than those predicted by 

the linear fit to the massive end of the relation, consistently with observations. This feature is not present if Supernova 

feedback is turned off. In the case of the baryonic Tully-Fisher relation, we also detect a weaker tendency for smaller 

systems to lie below the linear relation described by larger ones. This behaviour arises as a result of the more efficient 

action of Supernovae in the regulation of the star formation process and in the triggering of powerful galactic outflows 

in shallower potential wells which may heat up and/or expel part of the gas reservoir. 

Key words, galaxies: formation - galaxies: evolution - galaxies: structure 



1. Introduction 



The origin of the Tully-Fisher relation (jTullv fc Fisher I 
Il977f ) has been analysed by numerous observa- 
tional and theoretical works since it provides im- 
portant co nstraints for galaxy formation mod- 



els (e.g. Avila-Reese, Firmani. fc Hernandez I 



iMo. Mao, fc White I Il998t lAvila-Reese et al. I 



Originally defined as the relation between the luminosity 
and the rotation velocity for spiral galaxies, it is now 
accepted that this is a proxy of a more fundamental 
one between stellar mass and rotation velocity (e.g. 
Conselice et"an 120051: iFlores et al.1 120061: ICresci et al.1 
2009f ) or even between bary o nic m ass 
velocity (e.g. jBell fc de Jong l2001t_ 



and 



Gurovich et al. 
20101) of the" 



20041: 
form M 



Verheiien 



rotation 
200lt 



McGaugh I l2005t iPuech et al. I 



oc V^. Observational evi- 
dence suggests the existence of a break in the stellar 
Tully-Fisher relation (sTFR) at a rotation velocity of 
QOkms"^ (iMcGaueh et al. Il2000l ). Even more, recently 



iMcGaugh et al. I ( 2010f) also reported a break in the 
baryonic Tully-Fisher relation (bTFR) but at a lower 
velocity. 

Assuming that galaxies formed in virialized haloes and 
that the fraction of gas transformed into stars is inde- 
pendent of halo mass, a relation between stellar mass 
and circular v elocity with a slope of ct ~ 3 is di- 
rectly obtained ([White fc Frenk lllQQlHMo. Mao, fc White I 
[l998). This theoretical prediction is reproduced in cos- 



mological simulations when non effective mechanisms 
to regulate the gas c ooling and star formation ac- 
tivity are in troduced ("Tissera. Lam bas fc Abadil 119971 : 
ISteinmetz fc N avarro 1999). Therefore, changes in the 
slope of this relation to match observed reported values, 
which are on average close to a ~ 4, might require the ac- 
tion of other physical processes. Theoretical findings sug- 
gest that the triggering of powerful mass-loaded galactic 
winds by massive Supernova (SN) explosions would have 
an efficiency which depends on the virial halo so that 
large r effects are expected i n low mass systems (jLarson I 
I1974D . In fact. iDekel fc Silk I (|1986D estimated analytically 
that haloes with circular velocity smaller than r^ 100 km 
s~^ would be strongly affected by massive SN events. 
Nonetheless, the role of galactic outflows on the origin 
of the relations between the stellar and dynamical prop- 
erties of dwarf galaxies stil l constitutes an op en problem 
(jMac Low fc Ferrara I [1991 iTassis et al. I I2008D . In a cos- 
mological framework, only recently it has been possible to 
start simulating th e triggering of galactic outflows on more 
physical basis (e.g. lScannaoieco et al. I[2006t IStinson et al. I 
2007). In fact, cosmological simulations that include this 
process obtain ed sTFRs in more qualitative agreement wit h 
observations ([Governato et al. 1 120071: ICroft et al. 1 I2009D . 
There is evidence that the implementation of a suitable 
feedback model could account for a bend in the TFR as 
suggested by semi - analy ti cal works (e.g. Kang et al. 200j t 
Nagashimaet al I 120051: iDutton fc van den Boschi 12001 
Guo et alTI 20101) and some hydrodynamical simulations 
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(e.g. iTassis et al. 1120081) . However, a detailed analysis of 
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Table 1. Cosmological Hydrodynamical Simulations stud- 
ied in this paper. 



Name 


F/NF 


N, 


Mdark 

[lO^Me/i-i] 


[10 


A/gas 
«M0fe-^] 


S160 
S230 
S230NF 
S320 


F 
F 

NF 
F 


2 X 160^ 
2 X 230^ 
2 X 230^ 
2 X 320^ 


17.6 
5.93 
5.93 
2.20 




2.71 
0.91 
0.91 
0.34 



Col. 1: Name of the simulation. Col. 2: Model with/without 
SN feedback (F/NF). Col. 3; Initial number of particles in the 
simulation. Col. 4: Mass of dark matter particles. Col. 5; Initial 
mass of gas particles. 



the effects of SN feedback on the TFR using a physically 
motivated model was still missing. 

In this paper, we present results on the analysis of the 
sTFR and bTFR by using a SN feedback model specially 
suitable t o study the formation of galaxie s in cosmological 
scenarios (jScannapieco et aLll2005lT2006D since it does not 
include any scale-dependent parameter and yet, it is capa- 
ble to produce the expected anti-correlation between the 
strength of the galactic outflows and the p o tentia l well of 
the haloes, as shown bv IScannaoieco et al. I (|2006f ). 

This paper is organized as follows. Section [2] describes 
the numerical experiments. Section[3]outlines the results. In 
Section2]we discuss about our main flndings and in Section 
[5] we present the conclusions of this work. 



2. Numerical Experiments 

We performed cosmological simulations of a typical field re- 
gion of the Universe consistent with the concordance model. 
We adopt the following cosmological parameters: ft^ — 0.7, 
i7m = 0.3, fib ~ 0.04, a normalization of the power spec- 
trum of tJs — 0.9 and Hq — 100 h km s~^ Mpc"'^, with 
h = 0.7. Two of the simulations have been ran with 2 x 230'^ 
(S230) particles while a third and a fourth one used 2 x 320'^ 
(S320) and 2 x 160^ (S160) to check numerical effects. In ah 
cases, the simulated volumes correspond to a cubic box of a 
comoving 10 Mpc /i~^ side length. The particle masses for 
the dark and initial gas components are given in Table [1] 
Simulation S320 could only reach z ~ 2 due to high compu- 
tational costs. By analysing the sTFR and bTFR in S230, 
S320 and S160, we check our results for numerical effects. 

The simulations have been performed by using a version 
of GADGET-3, an update of GADGET-2 optimized for mas- 
sively parallel simulations of highly inhomoge neous systems 
(|Springel fc Hernguist I [20031 ISpringel Il2005[ ). This version 
of GADGET-3 includes the multiphase model for the in- 
terstellar medium ( ISM) and t he SN feedback scheme of 
IScannapieco et aU (|2005l I2006D . We run the same initial 
condition with (S230) and without SN feedback (S230NF) 
in order to quantify the effects of SN-induced galactic 
winds. 

The SN feedback scheme adopted in this work considers 
both Type II and Type la SNe for the chemical and energy 
production. We adopt 0.7 x 10^^ erg per SN event. The 
chemical e volution model used i n this code is the one de- 
veloped bvlMosconi et al. I (|200lt) and, later on, adapted by 
IScannapieco et al. I (120051 ) for GADGET-2. This model as- 
sumes the instantaneous recycling approximation (IRA) for 



SNifl. Lifetimes for the progenitors of SNIa are randomly 
selected in the range [10^,10^] yr. The chemi cal yields for 
SNII are given by Wooslev fc Weaverl (I1995D while those 
of SNI a correspond to the W7 model of [Thielemann et al. I 
(IT991 . 

The SN feedback model is grafted onto a multiphase 
model specially designed to improve the description of the 
ISM, allowing the coexistence of diffuse and dense gas 
phases. Within this new framework, the injection of en- 
ergy and heavy elements into different components of the 
ISM can be more realistically treated leading to a more ef- 
fective production of galactic winds. The ejected energy is 
distributed into and managed differently by the so-called 
cold (temperature T < Tc where Tc — 8 x 10^ K and 
density p > O.lpc where p^ is 7 x 10~^^g cm~^) and hot 
(otherwise) ISM components. The hot phase thermalizes 
instantaneously the energy while the cold phase builds up 
a reservoir until it accumulates enough energy to raise the 
entropy of the gas particle to match the value of its own 
surrounding hot environment. The properties of the hot 
and cold gas components of each gas particle are estimated 
locally without using any global properties of the systems. 
The fraction of SN energy distributed into the cold phase 
is given by the feedback parameter ec- The ejection and 
distribution of chemical elements is coupled to the energy 
procedure so that the fraction of metals pumped into the 
cold and hot phases are also given by Cc in these sim- 
ula tions. We adopt Cg = 0.5 which have been found by 
Sca nnapieco et alTI (|2008D to best reproduce a disc galaxy 



in simulations of a Milky- Way type halo. 

Virialized structures are selected from the general mass 
distribution by using a standard friends-of-friends tech- 
nique. The substructures residing within a given virialized 
ha lo are then indi v idual ized with the SUBFIND algorithm 
of ISpringel et al. I (|2001[ ) to build up galaxy catalogues at 
different redshifts. Note that, by definition, a substructure 
can correspond either to the central galaxy or a satellite 
system within a given halo. In order to characterize the 
simulated galaxies, we define the baryonic radius (i?bar) as 
the one that encloses 83 per cent of the baryons associated 
to each substructure. In the case of S160, S230 and S230NF, 
for all the linear regressions, we only consider those simu- 
lated galaxies defined by more than 2 x 10^ particles within 
^bar, which is equivalent to have stellar masses larger than 
IO^Mq h-i. In S230, Milky- Way type galaxies at z = 
are resolved with approximately 10^ total particles within 
-Rbar- A similar procedure was carried out for the high reso- 
lution simulation S320, where Milky- Way type systems are 
already resolved with more than 10^ particles within -Rbar 
at z = 2.0. In S320, the smallest simulated galaxies we 
considered for calculations have lO'' particles within i?bar- 

The properties of the simulated galaxies such as stel- 
lar mass, baryonic mass and rotation curves are estimated 
within i?bar- We define as disc systems those which have 
more than 75 per cent of their gas component on a ro- 
tationally supported disc structure by using the condition 
a/V < 1 to select them (where a and V are the velocity 
dispersion and the tangential velocity component, respec- 
tively). We acknowledge the fact that, at z = 0, most of our 



^ IScannapieco et al. I (2005) also implemented a version that 
relax the IRA assu ming the age-mass-m etallicity fitting polyno- 
mials estimated by (|Raiteri et al. |[l996l ). However, no significant 
changes were found in the results compared to the IRA. 
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galaxies contains large stellar bulges and thick stellar discs, 
being more consistent with early type spirals. Nevertheless, 
the gaseous disc components are very well-defined and trace 
remarkably well the potential well of their host haloes as 
can clearly be appreciated in Fig. [1] Rotation curves are 
estimated by using the tangential velocity of gas particles 
on the plane perpendicular to their total angular momen- 
tum. For these systems, tangential velocity constitutes a 
good representation of the potential well of the system. For 
the sake of simplicity, we use V^^^ = GM{r < R)/R (where 
M{r < R) is the enclosed total mass within R) estimated at 
i?bar as an indicator of the gravitational mass for systems 
in both S230 and S230NF. 

3. Results 

We calculated the sTFR for our simulated galaxies 
by the fitting of a relation of the form logM* — 
Q;log(V'cir/100 km s^ ) + Yioo- In the case of the local 
sTFR in S230, we found a slope of 3.68 ± 0.09 and 
^100 = 9.42 lb 0.26 in g e neral agreement with observations 
(|McGaugh et al. Il200a [Bell fc de Jo^l200lD . As can be 
seen from Fig. [2] (upper left panel), the residuals of the 
linear fit to the relation (small box) depart systematically 
from zero at around 100 km s~^, so that systems with 
lower circular velocities tend to have smaller stellar masses 
than the predictions obtained from the fitting. Interestingly, 
this trend is con s istent with the observations reported by 



Table 2. Results of the linear fits to the high-mass end 
(M* > IO^Mq/i-I) of the sTFR for S230 at z = 0. 



iMcGaugh et al. (l2000h. and from theoretical expectations 
(jLarson ill974t iDekel fc Silk Ill986[) . In order to analyse at 



which extent SN feedback might be responsible for this be- 
haviour, we compared these findings with the ones obtained 
for S230NF as shown in the upper right panel of Fig. [5J We 
can appreciate that SN feedback seems to be crucial to 
reproduce the observed features in the sTFR since when 
this mechanism is turned off, the simulated sTFR exhibits 
a linear behaviour with a slope of ^ 3 in agreement with 
theoretical predictions. We can also see that the wind-free 
run predicts larger stellar masses for systems at a given cir- 
cular velocity since it cannot regulate the transformation 
of gas into stars which is too efficient in hydrodynamical 
simulations. 

We also studied the bTFR in S230 and S230NF (Fig. [H 
lower panels) . The bTFR is obtained by adding all baryons 
within i?bar, regardless of its physical state. We can see that 
both models (with and without SN feedback) predict a lin- 
ear trend for the bTFR, at least, for the range of velocities 
covered by these simulations. At z = 0, the simulated bTFR 
in S230 has a slope of 3.23 ± 0.08 and Yiqq = 9.56 ± 0.23. 
These values are in general good agreement with obser- 
vational results reported for late and early type galaxies 
(|De Riicke et al. 2007: G urovich et al. 2010). However, we 
warn that the comparison with observations is tricky since 
we sum up the total amount of gas, while observers are lim- 
ited by the instrumental techniques which give only access 
to gas mass with certain physical properties. In the case 
of S230NF, the bTFR has a larger value of Fioo because 
these galaxies have been able to retain their baryons in the 
central regions since no SN-driven outflows could be trig- 
gered. We also note that S230NF yields similar fittings to 
the sTFR and bTFR, indicating that stars dominate the 
baryonic phase at this redshift. 

As we have already mentioned, the previous analysis 
has been done by using the value of Vdrc at i?bar as a kine- 



Vi 


Yloo 


a 


Vcirc at 0.5i?bar 


9.51 ±0.26 


3.62 ±0.09 


\4irc at -Rbar 


9.42 ± 0.26 


3.68 ± 0.09 


Kirc at 1.5i?bar 


9.43 ±0.25 


3.59 ± 0.08 


V max 


9.41 ±0.25 


3.67 ±0.08 



Col. 1: Kinematic indicator used for the calculations. Col. 2: 
Stellar mass corresponding to a velocity of 100 km s~^ according 
to the fitting. Col. 3: Slope derived from the linear regression. 

Table 3. Results of the linear fits to the high-mass end 
(Af, > IO^Mq/i-I) of the bTFR for S230 at z^O. 



V, 


Yloo 


a 


Kirc at 0.5i?bar 


9.64 ± 0.20 


3.19 ±0.07 


Vcirc at -Rbar 


9.56 ±0.23 


3.23 ±0.08 


Kirc at 1.5i?bar 


9.57 ±0.25 


3.14 ±0.08 


V max 


9.55 ±0.24 


3.22 ±0.08 



Col. 1: Kinematic indicator used for the calculations. Col. 2: 
Baryonic mass corresponding to a velocity of 100 km s~^ ac- 
cording to the fitting. Col. 3: Slope derived from the linear re- 
gression. 



matical indicator. In order to analyse the dependance of 
our results on the particular choice of the radius at which 
we measure Vdrc, we also compared the relations obtained 
by using Vcirc evaluated at 0.5i?bar and 1.5i?bar- We also 
estimated the maximum value of the rotation curve T^nax 
as can be seen from Fig. [3] The simulated sTFR and bTFR 
are not significantly affected by variations in the radius at 
which V^circ is estimated. The largest changes with respect 
to our previuos results (Fig. [2]) are obtained when using 
V^irc at 0.5i?bar as the velocity estimator. Nevertheless, all 
changes remain whithin a a. It is also worth noting that, 
in these simulations, Vcirc at i?bar constitutes a good proxy 
for the velocity V^ax, which is commonly employed in many 
observational works. Tables [5] and [3] summarize the param- 
eters corresponding to the linear fits to the high-mass end 
(M* > lO^M0ft.-i) of the TFRs shown in Fig. H We can 
appreciate that, in all cases, the values of a and Yioo agree 
within a a. Regarding the characteristic velocity where the 
sTFR bends, our findings suggest that it does not depend 
on the particular kinematical estimators tested in this pa- 
per (Fig. [21 small boxes). In the light of these results, here- 
after we will continue using Vcirc at i?bar as the kinematical 
indicator for our calculations. 

By comparing the sTFR and bTFR in S230, it is clear 
that adding the gas mass in the calculations contributes 
to restore the linearity of the TFR over a larger velocity 
range. These findings suggest that, at z '^ 0, the bend of 
the local sTFR might be partially caused by a decrease of 
the star formation rate in smaller systems as a consequence 
of the gas heating by SNe. To quantify at which extent 
SN feedback can generate a decrease in the star formation 
activity of simulated galaxies, we calculated the star for- 
mation efficiency (eSFR) defined as the ratio of the star 
formation rate of a given galaxy to the total amount of gas 
contained within it. In Fig. 21 we compare the mean eSFR 
as a function of V^jic for galaxies in S230 and S230NF at 
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Fig. 1. Projected gaseous mass distributions for a typical disc galaxy in a face-on (left panel) and an edge-on (middle 
panel) views and the corresponding Vrot (dots) and X4;irc (continuous line) for this system (right panel) in S230. For this 
galaxy, i?bar is ~ 13.7 kpc h~^. 




1.6 1.8 2.0 2.2 

log(Vc,,c) [km s-'] 
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a = 2.96±0.02 i 
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Fig. 2. The sTFR (upper panels) and bTFR (lower panels) in the S230 (left) and S230NF (right) runs, at z = . Grey 
symbols depict low mass systems (Af* < lO^Mo/i"^), while massive ones {M^, > IO^Mq/i"^) are coloured in black. 
Systems in S230 run, with well-defined gaseous discs are enclosed with squares. The black line corresponds to the linear 
fit. The values of the slope (a) and Fioo are indicated in the figure. The vertical lines indicates the velocity where 
the residuals (small boxes) of the sTFR in S230 depart systematically from zero. The thicker lines denote the range of 
velocities associated with stellar masses larger than IO^Mq/i"^. 
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Fig. 3. Mean sTFR (left panel) and bTFR (right panel), and the corresponding standard deviations for S230 at z = 0. 
Results obtained by using different kinematic indicators are compared: K;irc at 0.5i?bar (blue), -Rbar (black) and 1.5i?bar 
(red), and the maximum value of the rotation curve Vmax (green). The dotted black lines represent the fittings to the 
high mass-end (Af* > IO^Mq/i"^) of the simulated TFRs when using Vdrc at i?bar as kinematic indicator. The vertical 
line depicts the characteristic velocity where the sTFR bends. The small boxes show the residuals of each TFR with 
respect to the linear fit to its high-mass end. 
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Fig. 4. Mean star formation efficiency (eSFR) as a func- 
tion of Fcirc for galaxies in S230 (triangles) and in S230NF 
(diamonds) at z = 0. 



z = 0. As expected, the model which includes SN feedback 
predicts smaller eSFRs at a given Kirc, with slow-rotating 
systems exhibiting the most important variations, on aver- 
age. These systems still have an important fraction of gas 
but it is in the form of a diffuse warm environment which 
does not fulfil the conditions to form stars. 

To analyse the evolution of both the sTFR and bTFR as 
a function of redshift, we performed linear fits to the fast ro- 
tators using the characteristic velocity Vdrc — 100 km s^ 
as a reference value at each analysed redshift. As can be 
seen in Fig. [SJ the residuals of the sTFR show a system- 



atic departure from zero already from z « 3 which oc- 
curs at approximately a similar characteristic velocity of 
~ 100 km s~ . For this SN model, we detect an evolution 
in the yioo of the sTFR of « 0.44 dex since z ~ 3, while the 
slope of the fast rotators systems remains almost constant. 
In the case of the bTFR, Yiqo evolves by around ^ 0.30 dex 
between z = 3 and z = 0. In particular, as can be appre- 
ciated from Fig. [5j the simulated bTFR retains the linear 
relation over a larger velocity range. At z = 0, there is no 
clear signal for a change of the slope within the numeri- 
cal dispersion. This finding is not in disagreement wit h the 
recent observational results of iMcGaueh et alH (|2010[ ) who 
reported a break in the bTFR at ~ 20 km s~^, because 
our simulations do not resolve numerically such low veloc- 
ity systems. While larger statistics and higher numerical 
resolution is needed to test the behaviour of the bTFR for 
low velocity systems, more detailed observations are also 
required to robustly probe the existence of the bend in the 
bTFR. With regard to the S230NF, both the sTFR and the 
bTFR exhibit a linear trend since z = 3 with no significant 
changes in the slope but with an evolution of « 0.55 in 
Yioo, which is consistent with the higher eSFR measured in 
the galaxies in this run. Both TFRs show larger dispersion 
for low velocity rotators which are also more gas rich than 
systems at the high velocity end. At a given velocity, slow 
rotators can have different baryonic and stellar masses ev- 
idencing their different evolutionary paths. This behaviour 
can be also seen in S230NF albeit weaker. Hence, the dis- 
persion has at least two causes: the different evolutionary 
paths at a given circular velocity which regulate the trans- 
formation of gas into stars and the action of SN feedback. 

To assess the effects of numerical resolution on the de- 
termination of characteristic velocity Vc ~ 100 km s~^, we 
performed linear fits to the massive end (Af* > 10^MQh~^) 
of the sTFR and bTFR in S320, S230 and S160 at z = 2.0. 
In Fig. ini we compare the mean values and standard de- 
viations of the residuals associated to each simulation. We 
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Fig. 5. The evolution of the mean sTFR (upper panel) and 
bTFR (lower panel) as a function of redshift for S230 for 
z = 3, 2, 1 and (green, magenta, blue and red). The dotted 
black lines correspond to the fittings to the simulated TFRs 
for systems with Vdrc > 100 km s~ at z = 0. The verti- 
cal line indicates this characteristic velocity. In the small 
boxes we included the residuals estimated with respect to 
the corresponding linear regressions calculated at each red- 
shift. The thicker lines in the residuals indicate the veloci- 
ties associated with stellar masses larger that 10^Mq/i~^. 



also show with symbols the results for S320. As can be seen, 
the trend for the sTFR is clearly present in all simulations, 
showing a change in the slope at a similar velocity. With 
respect to the bTFR, it also exhibits a bend in the residuals 
at around ~ 100 km s~^ at z '^ 2 but the trend is weaker 
than that of the sTFR. The three simulations agree on the 
existence of the departure from linearity at approximately 
the same velocity. In the three runs, we also find larger dis- 
persions for low velocity rotators. The agreement between 



Fig. 6. Residuals corresponding to the linear fits to the 
massive end (Af* > IO^Mq/i"-^) of the sTFR (upper panel) 
and bTFR (lower panel) in S320 (symbols) at z = 2.0. Blue 
symbols depict low mass systems (M* < IO^Mq/i"-'^), while 
massive ones (M* > IO^Mq/i"-'^) are coloured in black. 
Systems with well-defined gaseous discs are enclosed with 
squares. The solid lines depict the mean values and stan- 
dard deviations associated to the residuals of the linear fits 
of S320 (black), S230 (red), S160 (green) and S230b (cyan). 

the three different resolution-level simulations suggests that 
these results are robust against numerical artefacts. 

Finally, since the SN feedback model uses Tc = 8 x 10^ 
K (equivalent to log Kir ~ 1-95, where Kir is the virial 
velocity of the halo) as a parameter to separate cold and hot 
gas-phases surrounding stellar populations, we evaluated at 
which extent this selection can infiuence the determination 
of the characteristic velocity where the sTFR bends. To 
assess this issue, we performed an extra simulation (S230b) 
with the same initial conditions and SN and star formation 
parameters of S230, except for Tc which was lowered to 
5 X 10* K (equivalent to log Kir ~ 1.84). Fig. JHlshows that 
for both TFRs the residuals of these simulations are very 
similar. Particularly, both simulations predict a bend at 
^ 100 km s~"^ independently of the adopted value of Tc 
indicating that this parameter has no significant effect of 
the determination of the bend. 

4. Discussion 

To further investigate our findings regarding the origin of 
the bend of the TFR and the role played by SN feedback, 
we studied the fraction fi — fini-^i/^b-^vir for each simu- 
lated galaxy at z = 0, where Mvir is the total mass within 
the virial radius and i denotes the stellar or the baryonic 
component within i?bar- Hence, the quantity /* (fb) rep- 
resents the ratio between the stellar (baryonic) mass in a 
given galaxy and the expected baryonic mass within the 
virial radius inferred from the universal baryonic fraction 
(r2b/^m) corresponding to the adopted cosmological pa- 
rameters. As can be appreciated from Fig. [71 /* is an in- 
creasing function of the circular velocity which, on its turn. 
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is a measure of the potential well of the systems. The small 
systems tend to have /* < 0.1 at z = 0, while for more mas- 
sive galaxies, /* reaches ~ 0.3. These findings are consistent 
with the fact that SN outflows are more efficient to heat up 
and/or expel the gas content from star forming regions in 
shallower potential wells, leading to a decrease in their stel- 
lar mass content and making them to lay below the sTFR 
determined by massive galaxies. At z — 2, the trends are 
similar to the local ones with faster rotators systems ex- 
hibiting stellar fractions of around 0.4 and most of the slow 
rotators not exceeding 0.2. Nonetheless, it is worth noting 
that, at a given circular velocity, the discrepancies between 
the simulated stellar components and the theoretical bary- 
onic masses are larger in the local Universe, evidencing the 
accumulated effects of SN feedback. We also note that high 
redshift galaxies show a higher slope for the relation be- 
tween f.^, and Vc indicating the more important action of 
SN winds on smaller systems at this epoch. 

When the gas mass within i?bar is incorporated in the 
calculations, we obtained the total baryonic fraction /b. At 
z = Q, /b shows a weaker dependence on circular velocities 
than /*. More massive galaxies exhibit /b ~ /*, indicating 
that these systems are gas poor. On the other hand, small 
galaxies exhibit baryonic fractions /b which are more than 
a half larger than /* . These correlations explain the trend 
to recover a single slope relation for the local TFR over a 
larger velocity range when considering the baryonic content 
of galaxies instead of their stellar mass. However, the fact 
that /b < 0.55 within i?bar for our whole sample indicates 
that galactic outflows are important over the whole simu- 
lated sample. The effects on the faster rotators are also in- 
direct since they grow by the accretion of smaller substruc- 
tures which are strongly affected by SN feedback as shown 
before. It is worth noting, that the /* and /b predicted by 
simu lations are in quahtativ e agreement with observations 
(e.g. iMcGaugh et al. |l2010f ). albeit with a weaker depen- 
dence on M*. At z = 2, on the other hand, we obtained 
a stronger correlation between /b and T4irc- Slow-rotating 
systems at this redshift show, on average, similar values 
of /b to local ones. Conversely, fast-rotators in the local 
Universe exhibit half of the baryonic fractions derived for 
systems of similar circular velocities at z = 2. This be- 
haviour is consistent with the stronger signal for a bend 
in the bTFR that we found at high redshifts and shows 
that SN feedback has acted on these systems expelling part 
of the baryonic mass outside the very central regions from 
z-2. 

In order to estimate if the expelled baryonic mass is 
within the virial radius or if SN outflows have been pow- 
erful enough to transport it into the intergalactic medium 
(IGM), we calculated the fraction /^''' of the total baryonic 
mass within the virial radius of simulated galaxies rela- 
tive to the one derived from the universal baryonic frac- 
tion (fib/^m) in the adopted cosmology. The results can 
be appreciated in the right panels of Fig. [T] By compar- 
ing these findings with the ones obtained for /b, it is clear 
that an important amount of the missing baryons within 
-Rbar can be found in the surrounding halo. However, given 
that Z^''' < 1, there is also a fraction of the gas which has 
been blown away as a consequence of very efficient galactic 
winds. These effects are more evident at lower redshifts and 
for smaller systems. 

We can also see that the correlation with virial mass is 
clearly defined at both redshifts albeit with high dispersion. 
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Fig. 8. Fraction of stellar mass residing in simulated galax- 
ies at z = relative to the dark component within the virial 
radius. For symbol codes see Fig. [5] We show for compar- 
i son r esult s from the semi-anal ytical works by iGuo et al. I 
(|20100 and lMoster et al. I ()2009[ ). as indicated in the figure. 
Note that the vertical axis has been plotted using logarith- 
mic scale. 



At z = 2, the fraction of absent baryons within i?vir ranges 
from 20 per cent to 70 per cent, with slow rotators having 
experienced the most important losses, as expected. In the 
case of local galaxies, the fraction of lost baryons is within 
the range 30 per cent and 80 per cent. 

Finally, in Fig. |8l we contrast our findings with 
the trends obtained f rom the semi - empir ical models of 
iGuo et al. I (|2010D and iMoster et ahl ()2009f ) which are de- 
rived from the observed luminosity function. This figure 
shows the fraction of stellar mass (i^*) contained within 
each simulated galaxy at z = relative to the dark com- 
ponent within the virial radius. We see that there is a very 
good agreement between our results and previous works. 
However, we note that our model still tends to overpredict 
the stellar mass fractions at the low-mass end of the rela- 
tion. 

The dynamical evolution of galaxies in our simulations 
is driven by the joint action of several physical processes 
such as the star formation mechanism, gas infall and galac- 
tic outflows. The star formation process is triggered mainly 
when the gas gets cold and dense . As a consequence of 
star formation, SN energy is released heating up the sur- 
rounding cold and dense ISM which, when appropriated 
(see section [5]), is promoted to the hot phase. This gener- 
ates a stronger decrease in the star formation activity in 
systems with shallower potential wells because it is easier 
for the cold and dense gas in small systems to match the en- 



^ Recall that in this model the cold phase is defined as the 
gas component with temperature T <Tc where Tc = 8 x 10* K 
and density p > O.lpc where pc is 7 x 10~^^g cm~^ (Section[21). 
Otherwise, the gas is classified as hot phase. 
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Fig. 7. Left and middle panels: Fraction of stellar (/*) and baryonic (/b) mass residing in simulated galaxies relative to 
the expected baryonic mass within the virial radius for S230. Right panel: fraction of the total baryonic mass within the 
virial radius relative to the expected one {f^") for S230. Results are shown at z = (upper panels) and z — 2 (lower 
panels). For symbol codes see Fig. [21 Note that the vertical axes have been plotted using logarithmic scale. 



tropy of the hot and/or diffuse phase residing within their 
small dark matter haloes. But, the cooling times (tcooi) of 
the promoted particles are too short compared to their dy- 
namical times (tdyn) to allow this hot phase to be stable 
for a long time 0. Hence, the gas which has just reached 
the hot phase due to SN heating might cool down again in 
short timescales, returning to the cold phase. Therefore, SN 
feedback leads to a self-regulated cycle of heating and cool- 
ing exerting an important influence on the regulation of the 
star formation process in low-mass galaxies. As the systems 
get larger, the hot phase is established at a higher temper- 
ature due to the joint effect of the continuous pumping of 
SN energy into the hot environment and the increase of the 
virial temperature of the dark haloes hosting them. Hence, 
as the cooling times get longer compared to the dynami- 
cal time scales, the hot gas is able to remain in this phase. 
While the hot environment is able to build up as galaxies 
get more massive, galactic winds would be more difficult to 
be triggered in these galaxies, since more SN energy would 
be required by the cold phase in order to match the en- 



In our model, the cooling and dynamical times are calcu- 
lated following standard definitions as tcooi = 3kTfimp/2Ap and 
tdyn = y37r/6Gp, respectively. 



tropy of its nearby hot environment and, therefore, to be 
promoted generating outflows. Meanwhile, the cold gas re- 
mains available for star formation. It is also true that as 
the systems grows, their gas reservoir decreases, also de- 
creasing their star formation activity and consequently, the 
source of SN energy. Therefore, the action of SN feedback 
is much less important at regulating the transformation of 
the remaining gas into stars in these large systems. 

In order to quantitatively probe that these behaviours 
are taking place in our systems, we estimated the relations 
between the dynamical and the cooling times of the cold 
gas that has joined the hot phase by the action of SN feed- 
back (i.e. promoted particles), along the evolution of our 
simulated volumes. The fraction of gas-phase non-available 
for efficient cooling within simulated galaxies (F^^) is de- 
fined as the gas mass w ithin Ryj^ which satisfies the con- 
dition tcooi/idyn > 0.1 (Navarro fc White IfToOah . In Fig. 
[9l we plot F™ as a function of Vcirc in S230 at z = 
(left panel) and 2 = 2 (right panel). As expected, F™ 
is an increasing function of Vcirc- At z = 0, systems with 
V^irc < 100 km s~^ has an important fraction of their gas- 
phase (> 90 per cent) subject to an efficient cooling. Only 
at Khc > 100 km s~-^, F™ starts to increase reaching 0.8 
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at Vcirc ~ 250 km s~^. The same trends are present at 
z = 2, albeit with a larger dispersion at the high- velocity 
end. Interestingly, our model predicts that the transition 
from efficient to inefficient gas cooling occurs at around the 
characteristic velocity Vc ~ 100 km s"'^ where the TFR 
bends. This velocity is also in agreement w ith the theoret- 
ical expectations from iDekel fc Silk I (|l986fl who concluded 
that only for systems with Vdr larger than ^ 100 km s~ , 
the cooling times of the diffuse and hot gas are at least one 
order of magnitude longer than the dynamical times and, 
therefore, the gas could remain hot. 

5. Conclusions 

Wc have studied the sTFR and bTFR by performing hydro- 
dynamical simulations in a cosmological scenario including 
the actio n of the combined multiphase and SN feedback 
model of IScannapieco et al. I (|2005l 120061 ). Our results sug- 
gest the existence of a change in the slope of the sTFR and 
bTFR produced by the action of SN feedback. At z = 0, 
the simulated sTFR is very tightl y defined and in genera l 
agreement with observations Ce.g. iMcGaugh et al. Il2000h . 
We detect a systematical change in the slope of the sTFR 
from approximately ~ 100 km s~ , at least since z = 3, 
so that lower velocity systems tend to have smaller stellar 
masses than those predicted by the linear regression of the 
fast rotators. The comparison with higher and lower reso- 
lution runs shows that this feature is not just a numerical 
artefact. Regarding the local bTFR, over the velocity range 
resolved by our simulations, it can be fitted by a single lin- 
ear model. At higher redshifts, we detect a weak trend for a 
change in the slope. Similar trends are found by using Knax 
as the kinematical estimator. 

Our results suggest an evolution of « 0.44 dex from 
z = 3 to z = for t he sTFR which is gl obal agreement 
with that reported bv ICresci etHTI (|2009D of 0.41 ± 0.11 
dex from z « 2.2 to z = 0. For the baryonic relation, our 
model predicts a slightly weaker evolution of sa 0.30 dex 
from z = 3 to z = 0. In the wind-free run, both the sTFR 
and the bTFR show the same level of evolution (« 0.55 
dex) as expected since baryons are efficiently transformed 
into stars as soon as they collapse within the potential well 
of a galaxy. 

The characteristic velocity of ^ 100 km s~ appears 
naturally as a consequence of the physically motivated 
multiphase and SN feedback treatment adopted in these 
simulations. This model is successful at establishing a self- 
regulated star formation process without using any global 
properties of the systems but resorting to the local estima- 
tions of the thermodynamical properties of the gas in the 
distinct components of the ISM. We have checked that this 
result is robust against variations of the SN input parame- 
ters. Nevertheless, we acknowledge the existence of several 
issues that need to be solved in the future such as the level 
of change in the slope which seems to be weaker in our 
simulations compared to observations. In this regard, our 
findings suggest that the study of the sTFR and bTFR at 
different cosmic times could help to constrain SN feedback 
models. 
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Fig. 9. Fraction of the gas component within i?vir which satisfies the condition tcooi/idyn > O.f in S230 at 2 = 0.0 (left 
panel) and z — 2.0 (right panel). For symbol codes see Fig. [21 



